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Abstract 

The pseudorapidity measurements at LHC, although in the central region only, allows to perform preliminary tests of 
the multiparticle production extrapolation formula inspired by the recent cosmic ray data analysis. Feynman scaling 
violation in the form proposed originally by Wdowczyk and Wolfendale in 70s has been applied to the Pierre Auger 
Observatory and the Hi-Res group measurements. The consistency of the Extensive Air Shower development and 
anisotropy data was found for smoothly rise of the scahng violation parameter. We have shown that the longitudinal 
momenta of produced particles determined inclusively as rapidity (pseudorapidity) distributions measured by LHC 
experiments follow the some universal high energy distribution scaled respectively. The high degree of Feynman 
scahng violation is confirmed. The decrease of the very high energy interaction inelasticity suggested by cosmic ray 
data analysis is found to be consistent with LHC measurements up to 7 TeV. 
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1. Introduction 

The inclusive description of minimum bias LHC events is not as spectacular as, e.g., Higgs hunting, but is essential 
for other very important scientific endeavours. One of them is the Ultra High-Energy Cosmic Ray (UHECR) problem 
and the answer to the question of an existence of Greizen-Zatsepin-Kuzmin (GZK) cut-off 1 1 ] . The origin and nature 
of cosmic rays is studied for almost exactly 100 years. The great experimental effort has been taken recently by 
two groups: the Pierre Auger Observatory [5| and the Hi-Res experiment f^. The progress is observed, but the 
answers are still not decisive. The cosmic rays of energies of about lO^^eV, if they are protons, should not reach us 
from cosmological distances. On the other hand anisotropy measurements show that they probably actually do. Our 
knowledge about the nature of UHECR is based on observation of giant Extensive Air Showers (EAS) - cascades of 
secondary particles created in the atmosphere when the single atomic nucleus (proton in a simplest case) enters from 
above. It is expected that the EAS initiated by protons and iron nuclei should differ. This difference is determined 
by the rate of energy dissipation. Thus it depends strongly on the distribution of secondaries produced in the forward 
direction and on the nature of primary particle: its atomic mass. The long-lasting discussions on the primary cosmic 
ray mass composition at the very end of the cosmic ray energy spectrum, in the so-called "ankle" region ( isiab > 
10'^ eV), could not be conclusive also because of the lack on the more exact knowledge of the very high energy 
interaction physics, what makes the importance of the high energy proton fragmentation even greater for cosmic ray 
physicist, astronomers and cosmologists. 

Searching for regularities and phenomenological description of the multiparticle production model is as old as 
the modeling in high-energy physics itself. Starting from simple Fermi thermodynamical model, to the first parton 
(quark) model propositions by Feynman, the model extrapolation to much higher, cosmic ray energies was one of the 
most important and most wanted model predictions. It is usually in the form of a kind of scaling. The idea of limited 
fragmentation |4] applied to the quark-jet hadronization led to introduction of the Feynman scaling variable of Xf 
and the universal fragmentation function f(xf, s) - fpixp) Jsj]. This brilliant idea works well for the first collider 
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experiments up to ^fs ~ 60 GeV. However, when applied to cosmic ray EAS development, it was questioned already 
at the "knee" energies of Eiat ~ 10'^ eV. The SPS (^fs ~ 200 - 900 GeV) experiments allow to quantify the scaling 
violation. The scale-breaking model of Wdowczyk and Wolfendale has been proposed to described the CR data at 
the beginning of '70 |0]. It is, in a sense, a generalization of the Feynman scaling idea introducing the one scaling 
violation parameter 

In Ref. |7] we have shown that the light composition suggested by the studies of the anisotropy and the average 
depth of the shower maximum (xmax) does not contradict other results, mainly the width of the jCmax distribution, only 
if one assume strong Feynman scaling violation. 

The rapidity (pseudorapidity) distributions were measured by LHC experiments: ALICE||8i], CMS|@, [loll and 
ATLAS 01 in (the last for p± > 0.5 GeV only) in the central rapidity region I?;! < 2.5 for c.m.s. energies of 900 
GeV, 2.3 TeV and 7 TeV Narrow range of a rapidity (pseudorapidity) at first sight does not allow to study important 
characteristics of very forward particle production. To study the fragmentation region new measurements, specially 
by much forward detectors (LHCf), are welcome. But, as it will be shown below, the existing data can be used to test 
the scaling violation picture found in UHECR physics domain. 



2. Rapidity distribution 

Rapidity distributions measured in LHC experiments cover the central region where the produced particles are dy- 
namically separated from the valence quarks of colliding hadrons. The central rapidity density p(0) - l/cr (dcr/dy)|^ ^ 
is the variable describing the particle production there. The original Feynman scaling preserves the value of the cen- 
tral rapidity density. The plateau in rapidity is characteristic feature of independent jet fragmentation model as well 
as statistical models with limited transverse momentum phase space. Unfortunately, it is known for long, that such 
simple picture does not work. 

The phenomenological fit of the p(0) rise made more than twenty years ago in Ref. [12] is still valid. The 900 
GeV LHC measurements match well SPS UA5 result. The systematic discrepancy seen by CMS detector does not 
change this general opinion. 

2.1. Feynman scaling 

Feynman scaling Jstl can be expressed introducing one universal function ff of the variable x - p\\/pmdx which 
describes the invariant momentum (longitudinal py) distribution of particles crated in the high-energy inelastic (and 
non single diiTractive) interaction 

E 1 rfV 

= f(x, p±, s) = /f(x, pJ (1) 



^/s/2 cr dxd^pj_ 

where yfs is the interaction c.m.s. energy, E, py and p± are energy, and longitudinal and transverse momenta of 
outgoing particles (pmax ~ V^/2). Change of variable from Feynman x to rapidity y gives 

1 c/V 
cr dy d^p± 



where ^(y) = -^p^^T'm^/( V^/2) sinh(y). Using an approximate relation -^p^^T'm^sinh(y) a; pj^ sinh(?7) and intro- 
ducing the very convenient variable: pseudorapidity 77 = - In tan(0/2) we have 

-;T^=/.fesmh(,),p.) . (3) 
(T dT] d^p± \ y/s I 

The integration over all p± is obvious with uncorrected p^ and py and the universality of the p± distribution 

1 d(r /2<pJ . \ 

- — ^ Ffl — — smh(?7) . (4) 



o- drj \ ^fs 

The factor (p^) is a constant related to the transverse momentum scale. 
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Figure 1: Pseudorapidity distributions shifted by Ay = In( VJ/m) for ISR, SPS and Tevatron measurements (a), and distributions measured by LHC 
experiments at energies from 900 GeV to 7 TeV compared with SPS ^fs = 546 GeV UA5 result (b). 



We are interested of the extremely forward part of the (pseud o)rapidity distribution - projectile fragmentation 
region. It is convenient to move the longitudinal momentum distribution to the anti-laboratory frame (rj — > rj') 
where the projectile is at rest prior to the collision. This is done shifting the c.m.s. (pseudo)rapidity distribution by 
Ay = In ( yfsim) 



sinh(77') = sinh(77 - Ay) = sinh^;; - ln(-\As/m)) 



g,-ln(VI/m)/2 ]!L ^ ^ sinh(7,) . (5) 

After such transformation the direct comparison of particle production at different values of interaction c.m.s. energy 
is possible 



1 dcr (2{pJ . \ 

sinh(77) ^ F,j{t]) . 



cr drj' 



(6) 



This form of Feynman scaling was tested e.g. in ref. 111211 and it is found that it is valid only very approximately. We 
can see this in Fig.[T^, where previous millennium data are plotted as a function of the anti-laboratory pseudorapidity. 
The recent data from CMS |9, 10] and ALICE |8] are shown in Fig.IIJ). 

It is known that Feynman scaling is violated at least by the continuous increase of the central rapidity density what 
is easily seen in Fig.[T] 

2.2. Feynman scaling violation 

The original Feynman scahng implies that the inelasticity of proton-proton interaction, defined as a fraction of 
incoming energy carried by newly created particle, is universal, the same for all interaction energies. The first obser- 
vations suggested an attractive value of 0.5. The rise of some characteristics of the interactions (like, e.g., average 

or central rapidity density we mentioned above) makes the assumption about the constancy of the inelasticity 
not quite well justified. Introducing the multiplicative factor proportional to the observed rise of the rapidity plateau 
to the right-hand side of Eq.© we can try to recover a form of scaling. Applying this procedure the simplicity of 
the original Feynman idea is lost and the next correction for the rise of the average transverse momentum could 
be introduced here as well. We have used in the present work the average transverse momentum rise of the form 
{pi_) - 0.413 - 0.017 ln(i) + .00143 In^(i) shown in Fig. 4 of Ref. 1,10.1 . The additional inelasticity control parameter 
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Figure 2: Pseudorapidity distributions siiifted and transformed respectively adjusting ap for ISR, SPS and Tevatron measurements (a), and distri- 
butions measured by LHC experiments at energies from 900 GeV to 7 TeV compared with SPS ^fs = 546 GeV UA5 result (b). 



is an index in a power law multiplicative factor. These two modifications lead according to Eq.© to only slightly 
more complicated scaling formula 

1 dcr IsV' I2{pj \ 

- — = — Ffl — — smh(77) . (7) 
o- dj] \so) \ ^/s I 

We have used the UA5 data measured at -</so - 546 GeV c.m.s. energy [12] as a datum. The very accurate mea- 
sured NSD pseudorapidity distribution have been used as a definition of the universal Ff function. We adjusted the Of 
parameter value to minimize the discrepancy between Eq.(|7} scaling prediction and the distributions of pseudorapidity 
measured at different energies: from ISR to 7 TeV of LHC. The results are given in Fig.|2l 

Values of Of increase from ~ 0.05 found for ISR 53 GeV to ~ 0.11 at LHC 7 TeV. The increase is statistically 
not very significant, at least for the overall inelasticity, what will be discussed later The accuracy of the data scaling 
according to Eq.© can be estimated with the help of statistical tests. The;^'^ values for the ISR and SPS are of about 
X^/NDF a: 40/20. The systematic uncertainties of the Tevatron and LHC results makes the /NDF smaller but the 
overall tendency seen in Fig. |2] suggests strongly that proposed modification of the Feynman scaling is not a right 
solution for the extrapolation of interaction properties to the very high interaction energies. 

2.3. Wdowczyk and Wolfendale scaling 

It was shown in Ref. 101 that the almost forty years old modification known as Wdowczyk and Wolfendale (WW) 
scaling |6] could be still satisfactory used to scale the interaction properties to the ultra high (> 10'^ eV) cosmic ray 
energies. 

The original idea of the WW scaling 

f(x,p_i_,s) ^ (s/soy fwwixis/soT, Pj_) (8) 

is an extension of the Feynman fragmentation formula of Eq. ([1]) (the limit for a = 0) with the possibility to get the 
'thermodynamical limit' of n ~ i'^'* with a = 0.25. 

The WW model in its version of mid '80 has been successfully used for the EAS studies around 'the knee'. Its 
extension introducing partial inelasticities (energy fraction carried by specific types of particles), and the transverse 
momentum rise with interaction energy dependencies, as discussed above, gave better description of the production 
of different kinds of secondaries. As a result of this improvements the first power-law factor index was released and 
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Figure 3: Wdowczyk and Wolfendale scaling with both parameters a and a' adjusted to each experimental data set. 



gave an extra model parameter. This more flexible formula was applied, e.g., in Ref. where the agreement of the 
WW model predictions and the UA5 measured rapidity distributions was shown. It should be mentioned that original 
Wdowczyk and Wolfendale model gave a complete description of the multiparticle production process to be used 
mainly in EAS studies, so it contains such details as partial inelasticities, transverse momenta, semiinclusive properties 
etc. The fit shown in Ref. 1 12] is the effective, average description of inclusive data of rapidity (pseudorapidity) only. 
In the present work we explore the WW scaling of the form 



1 do- 
er di] 









1 i^WW 



71)- sinh(77) — 



ff-l/2 ^ 

(9) 



where ) is the average transverse momentum at the datum interaction energy ( yfso - 546 GeV). 

We have adjusted first both a and a' parameters independently to get the best scaling performance. Results are 
given in Fig. [3] 

Obtained values of a and a' are shown in Fig. |4^. Horizontal lines show results from Ref. lll2f|( solid for a and 
dashed for a', respectively). The thick solid broken line is the result for a of our UHECR analysis JtI]. It is seen that 
the predictions from Ref. [7] and the LHC data are consistent. Although the large uncertainties, which are result of 
limited rapidity range as well as possible systematics, do not allow for any stronger conclusions. 

We can, however, use the UHECR data analysis predictions for the values of a and test if results of the fit, with 
such reduced free parameter space, remains in agreement with the WW scaling. It can be seen in Fig.|5] 

The data description is not much worst than the one presented in Fig. [3] The constancy of the a' suggested by 
WW original papers and seen in Fig. |4^, still holds as presented as in Fig. |4j5. 



3. Inelasticity 

In Ref. itzt] it is found quite unexpected high energy behaviour of interaction inelasticity coeflicient. It was obtained 
as a result of the experimental suggestion that the composition of the UHECR is quite light, contains a significant 
proton fraction. The WW model with the strong Feynman scaling violation leads to continuous decrease of the energy 
fraction released to the secondaries produced in very high energy interactions. Eq.(|9) gives the inelasticity energy 
dependence 

/ \(a'-a) 

K(s) = Ko 1-1 , (10) 
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Figure 4: W&W scaling pai'ameters predictions for a (solid symbols and solid lines) and for a' (open symbols and dashed line) adjusted to the data 
(a), and values of a taken from the UHECR analysis (J and only a' used as a free parameter of the fit (b). 
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Figure 5: Wdowczyk and Wolfendale scaling results with a set to the UHECR analysis data and a' adjusted to each experimental data set shown 
as in the Fig. [3] 
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Figure 6: Inelasticity calculated with WW scaling assumption (filled symbols 
inspired a (Fig.|5)- 



■ circles for both a and a' adjusted (Fig[3) and squares for UHECR 



while for the modified Feynman scahng formula Eq.dT) it is 



K(s) = Ko . (11) 

In the Fig.|6]we have shown results of our analysis. Open symbols show the fast rise of the inelasticity for modified 
Feynman scaling formula. Even if the af follow the lower energy, smaller value, in the UHECR domain the saturation 
is expected. Filled symbols were obtained for WW scaling. The solid line gives the predictions from Ref. [7] obtained 



using UHECR data. The dashed line is the fit from Ref. 112|] of the WW scaling parameters to SPS data. The value of 
0.5 is also shown. 

The open symbols are for the modified Feynman scaling with af parameter. Solid line shows the UHECR data 



analysis prediction from Ref. [7J. Dashed line is the inelasticity fit from Ref. Ill2ll . The 'canonical' value 0.5 is shown 
by short dashed line. 



4. Summary 

We have shown that the minimum bias pseudorapidity distributions measured by LHC experiments can be very 
well described with the scale-breaking Wdowczyk and Wolfendale formula. 

The scaling violation observed for the energies up to SPS yfs = 900 GeV and 1800 GeV in Tevatron was uphold 
recently in the analysis of new UHECR data. 

The phenomenological model of Wdowczyk and Wolfendale introduces two model parameters. The value of one 
of them: a, was originally found to be equal to 0.13 using interpolation of the Xf - p\\/pmdx distributions between 
Vs ~ 10 GeV and ISR energies. Later interpolations including SPS data gave the value of 0.18 and finally the effective 
value of 0.25 was found in Ref. |12]. The increase of the central rapidity density reported also in Ref. [12] suggests 
a - 2 X 0.105 = 0.21. This value gives the Extensive Air Showers development maximum position x^ax for proton 
initiated showers not far from measured llH [sl as it is shown in Ref. 0]. 
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The UHECR data suggests further smooth rise of the scale-breaking parameter. The first measurements at LHC up 
to 7 TeV c.m.s. energy agree with the trend observed at lower energies and seems to smoothly bridge accelerator re- 
sults and these on very high energy interaction of cosmic ray protons. The limited range of measured pseudorapidities 
does not allow for a stronger statement. The more forward particle production data is highly welcome. 

The rising inelasticity for (modified) Feynman scaling is obviously in contrary to the Wdowczyk and Wolfendale 
scaling and cosmic ray data. Comparing the pseudorapidity distributions in Figs.jSj) andjS}' we can say that the LHC 
pseudorapidity data analysis favours the second possibility. 
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